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Abstract 


We construct a relationship between integral and differential representation of second- 
order Jordan chains. Conditions to obtain regular potentials through the confluent super- 
symmetry algorithm when working with the differential representation are obtained using 
this relationship. Furthermore, it is used to find normalization constants of wave functions 
of quantum systems that feature energy-dependent potentials. Additionally, this relation¬ 
ship is used to express certain integrals involving functions that are solution of Schrodinger 
equations through derivatives. 
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1 Introduction 

A Jordan chain for a linear operator is a sequence of vectors that satisfies a particular system 
of equations. In the most elementary case, the linear operator is a square matrix and the 
corresponding Jordan chain is a sequence of generalized eigenvectors that are needed for the 
construction of Jordan’s normal form [I]. If the linear operator associated with the Jordan chain 
is a Hamiltonian, then the corresponding system of equations appears in the confluent version 
of the supersymmetry formalism (SUSY). The general purpose of SUSY is the construction 
of solvable quantum models (SUSY partners) and the manipulation of their spectra (spectral 
design). There is a large amount of literature on the topic, a selection of which can be found in 
the references of the reviews [213] or the more recent work [3] . Within the formalism of SUSY, the 
standard algorithm and the confluent algorithm are to be distinguished. Being well understood, 
the vast majority of applications uses the standard algorithm, while its confluent counterpart is 
much less known. An introduction to the algorithm and discussion of mathematical properties 
can be found in [M3 whereas some applications in EMI and references therein. In the 
confluent SUSY algorithm, partner models are essentially determined through Jordan chains. 
The system of equations associated with a Jordan chain is solved by transformation functions 
that are used to construct the SUSY partner to a given model. The solution of such a system 
of equations can be represented in integral form m or by means of derivatives E3- While the 
process of obtaining the integral form is more straightforward [lOl 18]. a closed-form evaluation 
of the integrals is often not possible. Since the differential approach is much more feasible in 


calculations 1191120] , it is of particular interest to find the mathematical relationship between 
the two approaches. In the present work, we determined the latter relationship and presented 
several applications regarding the confluent SUSY algorithm and the resolution of integrals to 
find normalization constants for energy dependent quantum mechanical systems. 

The paper is organized as follows. We introduce and discuss Jordan chains of second order 
in section 2. We present in section 3 two main application of the results obtained: the confluent 
SUSY algorithm and energy dependent potentials, showing that some integrals appearing both 
situations can be expressed through derivatives. Finally, our conclusions are presented in the 
last section. 

2 Jordan chains of second order 

Let D C M be an open interval. The following system of differential equations is referred to as 
Jordan chain of second order 

u xx (x, A) + [A - V{x, A)] u(x, A) = 0, (1) 

v xx (x, A) + [A- U(x,A)] v(x,X) = [V\(x, A) - l]u(x,A). (x,X) € (U,M). (2) 

Here, the indices denote partial differentiation. We assume that the function V and its partial 
derivative V\ are continuous. Furthermore, the solutions u and v are required to be three times 
continuously differentiable. Since equation ([!]) resembles the quantum-mechanical Schrodinger 
equation, the function V and the variable A will be called potential and energy, respectively. In 
the particular case that the potential is independent of the energy, we have V\ = 0, such that 
the Jordan chain simplifies to 

u xx {x, A) + [A - V{x)\ u(x, A) = 0, (3) 

v xx (x,X) + [A-U(x)] v(x,X) = -u(x, A), (x, A) G (D, M). (4) 

These equations play an important role in the confluent algorithm of the quantum-mechanical 
SUSY formalism, as will be discussed in detail below. Let us now return to our more general 
system ©, m , where we are interested in the construction of solutions. Our starting point is 
the assumption that we know a solution u of the first equation © . From here we will distinguish 
two options for determining a solution v to the second equation ©. 

2.1 Integral representation 

Our first option is the variation-of-constants formula, which requires two linearly independent 
solutions u\ and u 2 of © in order to become applicable. For the sake of simplicity, let us make 
the following choices for these two functions 


ui(x,X) = u(x, A) 


u 2 (x,A) 



(5) 


Note that the integral form of u 2 is obtained through reduction of order [21]. Observe further 
that due to the choice ©, the Wronskian W U1>U2 of u\ and u 2 is equal to one. The variation-of- 
constants formula then gives the following particular solution vyc to the second equation © of 
our Jordan chain. 


v V c(x,\) = 

X X 

= ui(x,X) j u(t,X) u 2 (t,X) [1 — V\(t, A)] dt — u 2 (x, A) J u(t, A) ui(t, A) [1 — V\(t,X)]dt, 



recall that the function u appears in the nonhomogeneous term on the right side of ([2]) . We 
can replace the latter function, as well as u 2 , by their definitions shown in ([5]). We obtain the 
following expression 


vvc{x, A) 



u 2 (t,X)[l -V x {t,X)]dt 



I u 2 (t, A) [1 - V\(t, A)] dt 


After applying integration by parts to the first term on the right side, we arrive at the following 
integral representation of our solution vyc 


v vc (x,\) 


—u(x, A) 


X 



t 

J u 2 {s, A) [1 - V\(s, A)] ds 


u 2 (t, A) 


dt, 


( 6 ) 


The principal disadvantage of this representation is its limited applicability. Inspection of ([6]) 
shows that the double integral can only be evaluated in closed form for relatively simple functions 
u and V. For a solution u of the second-order equation (JT]) that is typically given through special 
functions, it is in general not possible to evaluate <[6|) . 


2.2 Differential representation 

We will now derive a representation of the solution v to (J2]) that does not contain any integral. 
To this end, we substitute the derivative v = u x into the left side of CD). This gives 

(u x )xx(x,X) + [A - V(x, A)]u A (a:, A) = 

= (' u\) xx (x, A) + {[A - V(x, A)] u(x, A)} a - [A - V(x, X)] x u(x, X) 

= {u xx {x , A) + [A - V(x, A)] u(x, A)} a - [A - V(x, A)] A u(x, A) 

= - [A - V(x, X)] x u(x, A). (7) 

Note that in the last step we made use of ©• Next, we observe that (ED is the same as our 
initial equation (f2j) . Consequently, we have shown that 

v DF (x, A) = u x (x, A), (8) 

is a particular solution of Q, where the index is an abreviation for differential formula. Since 
© is free of integrals, it is generally much easier to evaluate than ©, and as such more suitable 
for applications. 

While both the variation-of-constants scheme and the differential formula provide meth¬ 
ods for solving (f2|) , the respective particular solutions vyc and vdf are in general different 
from each other. It is therefore desirable to know the precise relation between these functions, 
such that one function can be expressed through its counterpart. In order to find this relation, 
let u\ and U 2 be linearly independent solutions of ED, the Wronskian of which is equal to one. 
Recall that such a choice can be made through the settings ©. Since both vyc and vdf are 
particular solutions of the nonhomogeneous equation (f2l) . their difference must be a solution 
to the homogeneous equation ED> which can be given as a linear combination of u\ and U 2 - 
Equivalently, there are two parameters d\ and cfe, such that the equation 


di(X) ui(x,X) + d 2 (X) u 2 (x,X) = v DF (x,X) - v vc (x,X), 


(9) 













is fulfilled. The right side of this equation is completely determined by the choice of u\ and u 2 . 
Therefore, the remaining task is to calculate the functions d\ and d 2 . Since we need a second 
equation in order to find these two constants, we take the partial derivative of ([9]) with respect 
to x, yielding 

di(A) (ui) x (x,\) + d 2 (X) (u 2 ) x (x,X) = (v DF ) x (x, A) - (v VC )x(x, A). (10) 

Hence, we now have an algebraic system of equations ©, cnD that is linear in d\ and d 2 . 
Application of Cramer’s rule and taking into account that W U1>U2 = 1, we arrive at the Wronskian 
representation of the solutions 

dl(A) ^^VDF— v VCi u 2 (^O: ; <^2 (A) W U \,VDF—VVC (^o. A). (U) 

Note that the constant xo € D can be chosen arbitrarily, since neither d\ nor d 2 depend on the 
actual value of the variable x. In summary, if the functions d\, d 2 are chosen as in (Hip . then 
the particular solutions vyc an d vdf of © are related to each other by means of the identity 
(|9j). Recall that the explicit form of the particular solutions can be found in ([6]) and ©. 


3 Applications 

The Jordan chain and the differential representation of its solution has several interesting appli¬ 
cations, two of which we will now introduce. The first application concerns the normalizability 
of solutions associated with Schrodinger equations for energy-dependent potentials. The second 
part of this section is devoted to the confluent SUSY formalism and several examples. 


3.1 The confluent SUSY algorithm 

The second application of our results on second-order Jordan chains concerns the confluent SUSY 
algorithm. For sake of completeness we present here a brief introduction to that algorithm, for 
a detailed description see [1] and references therein. 

We start out by considering the following stationary Schrodinger equation 

ipxx(x) + [e- V(x)\i/j(x) = 0, (12) 

where e is a real constant and the potential V is a real function. Suppose further that we know 
two auxiliary functions u, v satisfying the Jordan chain 


u xx (x, A) + [A - V(x)] u(x, A) = 0, (13) 

v xx (x, A) + [A- V(x)]v(x, A) = -u(x,X), (14) 


for a real constant A. Observe that these equations are a special case of our general system ©, 
m, where the potential does not depend on the variable A. Then the function 


4>{x,X) 


Wu,v,^ (xj A) 
1^11,11 (x, A) 


u 2 (x, A) 
}^u,v(xi A) 


1px(x) + 


u(x, A) u x (x, A) 


^(x), 


(15) 


when e / A, or in the opposite situation 


<Xx, A) 


u(x, A) 
W u , v (x, A)’ 


fulfills the equation 


(16) 


e — V(x, A) 


4>xx (x, A) Y 


4>{x, A) = 0, 


(17) 










where the potential V is given by the expression 

V(x, A) = V(x) - 2 {log [W u , v (x, A)]} M . (18) 

This transformation is known as confluent SUSY transformation and the potentials V and V are 
called SUSY partners. Observe that the function (j> in (1151) depends on e. Since we are mainly 
interested in the dependence on A, we do not include e as an argument of q A 

Now, as can be seen from (fT5l) . (fl6l) and (fl8l) the Wronskian W u>v plays a fundamental role to 
perform the transformation. Furthermore, to obtain regular potentials V we should avoid zeros 
in this Wronskian. If the transformation function v is obtained by the variation-of-constants 
method (see (j6])), then W u _ v can be expressed as 

W UtVvc (x,\)= uj 0 -[ u 2 (t,\)dt, (19) 

Jx 0 

where xq is a point in the domain of V and ojq is an arbitrary constant. The conditions for 
obtaining a regular potential V from the confluent SUSY transformation are known and can 
be presented as follows: Let the domain of V be the interval (x£,x r ), then we have two non¬ 
excluding options: 

• We use a solution of (1131) as transformation function such that u(xi) = 0. In this case the 
integral 

PXO 

h = I u 2 (y, \)dy < oo, (20) 

and as a consequence if coo E (— 00 , — Jg], we avoid zeros in the Wronskian. 

• The transformation function u satisfies u(x r ) = 0. Thus, the integral 

PXr 

I r = U 2 (y, \)dy < 00 , (21) 

Jx 0 

and if coq E [/ r ,oo), we guarantee that W U}VVC never vanishes. 

When the solution v of (|14l) is a superposition of the general solution of the corresponding 
homogeneous equation and the particular solution obtained by differentiating u with respect to 
the parameter A (see (j8|l) the regularity conditions have not been fully explored. Let us express 
the Wronskian by 


W UjV (x,A) = K + W u>ux (x, A), 


( 22 ) 


where K is an arbitrary constant. A derivative with respect to the position, and using the 
Jordan chain (1131) (1141) to simplify second order derivatives, leads to 


dW u , v (: r,A) 
dx 


-u 2 (x, A), 


(23) 


i.e., for a fixed A, W UyV is a non increasing monotone function of x. Thus, two conditions can 
lead to a regular potential: 


• Non diverging transformation function u(x, A) when x —>• xi and 


K < - W U!Ux (x£, A). 


(24) 


• Non diverging transformation function u(x, A) when x —>• x r and 


K > -W U}UX (x r ,\). 


(25) 



Note that when the transformation function fulfills both boundary conditions, the domain 
of K is (— 00 , -W UtUx (x£, A)] U [—W UjUx (x r , A), 00 ). 

Furthermore, integration of (1231) leads to the interesting result 

f u 2 (t, X)dt = W U:UX (xo, A) - W UtUX (x, A), (26) 

Jx 0 

where xq € (x£,x r ). This last equation shows how to integrate the square of a function that 
is solution of a Schrodinger equation using a Wronskian and vice versa. This identity is useful 
when we are interested, for example, in finding normalization constants or probabilities in an 
interval. 

In the following two examples we illustrate first with the infinite well quantum system how 
to obtain the normalization constant of its eigenfunctions using (1261) and then how to obtain its 
SUSY partners using vdf and the regularity conditions presented in (1241) and (051) . Since SUSY 
partners of this potential have been obtained in [22] using the representation of the Wronskian 
as in (usd, this example provides a good comparison of the two methods. In the second example 
we build new exactly solvable potentials departing from the radial oscillator system and study 
the regularity conditions for obtaining these potentials. 


3.1.1 Particle in a box 

Consider a particle included between two impenetrable potential walls at x = 0 and x = 1. The 
Schrodinger equation of the system is 

i/)xx(x) + e i(>(x) = 0, x € (0,1), (27) 

and the boundary conditions are 


lim ip(x, e) = lim i/j(x, e) = 0. 

x —>0 x —>1 

If e = k 2 the the general solution of (1271) is given by 

4>(x) = Alsin ( kx ) + B cos ( kx ) 


(28) 


(29) 


where A and B are arbitrary constants. To satisfy the boundary conditions (|28l) we need to set 
B = 0 and k = nir, where n is a natural number, then the eigenfunctions and eigenvalues are 

V’(x) = Alsin (mrx ), e = n 2 7r 2 , (30) 


the normalization constant A of which we can obtain by means of (1261) . In order to avoid 
confusion regarding the notation, we extend the solution (1291) to a function of two variables 
ij} = ^(x,e), where it is understood that \k\ = y/e. The derivative with respect to e can be 
obtained using the chain rule 




dk d-il’ 
de dk 


—— sin (kx) 
2k v ' 


(31) 


note that we did not include the variable of differentiation e as an argument of We have 




A 2 sin (2 kx) 
~2~ 2k 



finally, the normalization condition leads us to 


1 = 



Ww € (0,e) l,e) — 


(32) 


A 2 
2 


=> A = V 2 . 


(33) 







Figure 1: On the left, potential V using A = 47 t 2 and K = 0.555. On the right its first three 
eigenfunctions. 


To generate a exactly solvable potential through the confluent SUSY formalism using the 
differential representation of the Wronskian (|22D set A = m 2 ir 2 , where m is a natural number, 
then the transformation functions u and v satisfying the Jordan chain (1 1311 - (114(1 for V(x) = 0. 
Let us consider 


u(x, A) 
v(x, A) 



(34) 

(35) 


where K is an arbitrary constant. Note that the first term of v is a particular solution of (11411 
obtained as vdf = u\, while the second term is solution of the corresponding homogeneous 
equation. Now, to obtain regular potentials we need to evaluate the Wronskian W UiUx in x = xe 
and x = x r : 


W UjUx (0, A) = 0, W UiUx (l, A) = (36) 

According to the conditions (f2~il) and (1251) . if the constant K € (—oo,0] U [1/2, oo) then the 
produced potential is regular. The potential generated is given by (1221) and can be expressed as 


V(x,X) 


167 x 2 m 2 [1 + rmr(2K — x) sin(2m7rx) — cos(2m7rx)] 
[2mir(2K — x) + sin(27rmx)] 2 


(37) 


recall that m = s/X/tt. Figure Q] shows a potential V (see (1371) 1 where transformation functions 
(|34l) - (l35j) with a parameter A = 47 t 2 were used. On the right, it can be seen its first three 
eigenfunctions (continuous, dotted and dashed curves respectively), the first and third were 
generated with the rule (1151) while the second with (|16l) since A = 47 t 2 is also and eigenvalue of 
the original system, see (1301) . 

Moreover, any integral of the form f u 2 dx, where u is solution of ()13l) . can be obtained with 
(|26l) even though the boundary condition (|28l) are not satisfied. Let us consider 

1 

=-cos 

\/A 

then the following integrals can be obtained with (1261) : 



u\(x, A) = sin ( v Ax 


u 2 (x, A) 



- W uu{ui)x (x,\) = | - ^=sin (2VAx) , (39) 

^u2i(u2)a — ^u2,(u2)a Ty qy3/2 ^2\/Ax^ . (40) 














Furthermore, using the relations Q and (|11D doubles integrals of the form f (f u 2 dx)/u 2 dx can 
be found. The transformation function v using the variation of constant formula ([6]) in this 
example is given by 


vvc(x, A) 


r x r 

• x) j 

J Xn \_J X 


—u\{x 


u\(s, X)ds 


X 0 IJXQ 


ul(t, A) 


dt. 


(41) 


notice that the two integration constants associated with the integrals have been fixed by means 
of the integration limits, where xq is an arbitrary number in the domain D of our Jordan chain. 
Now, the function vyc satisfies the conditions 


vvc(x o, A) = (v V c)x(xq, A) = 0. (42) 

Keeping this in mind, the formulas in (11111 reduce to 

d\ (A) ^^V£)F,U2 (*^0) A); ^2 (A) ^^U1,VDF (*^0, A) • (43) 

Let us point out that the functions in (1431) do not depend on any integrals. Now, solving ([9]) for 
the double integral in it, 


rx r rt 

/ / Ui(s,X)ds 

Jx 0 Ux 0 


1 , 1 
-dt = — - 


u 2 (t, A) u±(x,X) 


[(rti) A (x, A) - di(X)ui - d 2 (X)u 2 (x, A)], (44) 


recall that vdf = ( ui)\■ Substitution of (1381) in (1U1) leads to 


rx I" rt 

/ / uf(s,X)ds 

j XQ VJ X 0 


COS' 

dt =- 


(\/Axo) 


«?(*, A) 


2A 


x — Xq 
2\/A 


+ 


sin (2VXx 0 


4A 


cot ( VXx ) . (45) 


In this simple case, integrals in (139|) . (fJUl) and (l45j) can be verified by direct integration. 


3.1.2 The radial oscillator system 

The Schodinger equation for the radial oscillator potential is given by 

’ + 1 ) 


l/>xx(x) + 


e — x — 


tp(x) =0, x € (0,oo), 


where the potential V can be identified as 


V{x) = x 2 + 


(46) 


(47) 


and we will consider in this example l as a natural number. The boundary conditions of this 
quantum problem are 


lim ijj(x) = lira i/j(x) = 0. (48) 

x —>0 x—>oo 

The general solution of (|46D without considering the boundary condition is 

ip(x) = x £+1 exp^-y^ | A 1 F 1 (^(21 + 3-e);£ + ^-x 2 ^j 

+ Bx-W + V \F\ Q(-2£+l-e);-£+i;x 2 H, (49) 

where i F\ is the confluent hypergeometric function [23], and A, B are arbitrary constants. Since 
the second term of '0 as given in (1491) is divergent at the origin we need to set B = 0 in order 





















to satisfy the first boundary condition. The second condition is fulfilled only when the first 
argument of the hyper geometric function is a negative integer. Thus, the eigenfunctions and 
eigenvalues of (1461) are 

ip{x) = A x e+1 exp | 

To obtain the confluent SUSY partners <118|) of the potential (I47|) using the differential formula 
(1221) we need a solution of the Jordan chain m (ED- We will use the solution u of the first 
equation in the Jordan chain as given in (1491) with A = 1, B = 0 and A different from any of the 
eigenvalues: 

u(x, A) = x i+1 exp i^i ^(2f + 3-A);f+^;s 2 j . (51) 

A particular solution of the second equation in the Jordan chain can be constructed as u\. In 
order to obtain this derivative, the formula of the derivatives with respect the parameters of 
the confluent hyper geometric function can be found in [25]. Then the second transformation 
function can be expressed as: 


\F\ ( -n\£ + -;x' 


e = 4n + 2£ + 3, n = 1,2,3,..., (50) 


v(x,X) = 


x £+1 

exp 

K 



x 

~Y 


oo 

E - 4 

m— 0 


(i(2£ + 3- A)) m x 2m 


V- 


S') 

2 


m\ 


m— 1 

, E 


Zn P + Z + 3 — A) 


p=0 


2 £ + 1 


exp (-y ) 1*1 ( \ (~ 2 ^+ 1 - A);-^ + 


(52) 


where (-) m is the Pochhammer symbol [ 23 ] and K is an arbitrary constant. Note that the first 
term in (1521) is the particular solution while the second is solution of the homogeneous equation. 
To simplify some expressions let us define an auxiliary function h = h(x, A) as: 

h(x,X) = ^(2l + 3) A ^ +2 6XP (~y) lFl (l^ + 7_ A ) ; ^ + l ;a;2 ) ’ 

Notice that the expression of the Wronskian W u>ux = uh\ — u\h. The expression of the SUSY 
partners V is in (]18[) and for the radial oscillator can be expressed as 

y, _ x 2 + + 1) 4tt(x, A )u x (x, A) [K + u(x, X)h\(x, A) - u\(x, X)h(x, A)] + 2u 4 (z, A) 

x 2 [K + u{x , A )h\(x, A) — u\(x, X)h(x, A)] 2 

According to the condition (1231) . and since JU UjUA (0, A) = 0, to obtain regular potentials the 
value of the constant K has to be a non positive real number. On the left of figure [2] it can 
be seen a potential V, see ([531) . and as a reference the radial oscillator potential (1371) (dotted 
curve). The parameters used are £ = 1, A = 8 and K = —0.01. On the right we plotted its 
first three eigenfunctions, the first and third (continuous and dashed curves, respectively) were 
obtained with the rule (fT51) while the second (dotted curve) with (fl6l) . 

Before we conclude this section, let us comment on one more application that arises within 
the present example of the radial oscillator system. As we observe in (1491) . the general solution 
of equation (1461) is a linear combination of two solutions that we will call u± and u^- Using the 
abbreviations 

a = - (2^ + 3 — A), b = £+~, c = - (— 2£ + 1 — A), d = —£ + ~, (55) 

the particular solutions in (1491) can be written as 

ui(x, A) = a3 +1 exp^-y^ (a; 5; x 2 ) , (56) 

u 2 (x,X) = - 2 / +1 x~ £ exp i-Fi (c; d; x 2 ) . 


(57) 













Figure 2: On the left, a potential V using a transformation function u with the parameters 
1 = 1, A = 8 and K = —0.01 (continuous curve) and as a reference the potential V (dotted 
curve). On the right its first three eigenfunctions. 


Since u\ and U 2 are solutions of equation (l46j) for e = A, the identity ([26]) can be used to find 
several integrals involving the latter functions. We obtain 


J u\{t,\)dt = exp (-z 2 ) <( iF) (a; 6; * 2 ) £ 


(a + l) n X 2 


E 


1 


(b+l) n n\ ^ p + a + 1 

n=0 x 7 p=0 


t? ! i u i i l^i (a;b-,x 2 ) ( a )n % 2n 1 

+ i-Fi (a + l;6 + l;x ) ---E m E 


(&)„ n! ^p + a 

n=0 v p=0 


, (58) 


and 


roo 

/ u\(t,X)dt 

J X 


C X 


l-2t 


2(2f+ Yfi exp W 2 ) ] ^ (« * O E ( ( 7 +VJ E 


p =0 


T i F\ (c + 1; d + 1; x 2 ^ 


\F\ ( c;d-,x 2 ) _ (c) n g 2n yl 1 

c h p+c 


• (59) 


Furthermore, since u = ui + U 2 is also solution of the same differential equation and due to the 
linearity of the Wronskian we can obtain the integral of the the product u\U 2 - This results in 

























the comparatively long expression 


px 

/ ui(t, X)u 2 {t, X)dt 
Jx 0 


8 ( 2 £ + 1 ) 

1 F 1 (c; d;x 2 ) Y 


(2£ + 1) exp (—x 2 ) 
‘ (a)«x 2nn " 1 


°° ( \ 2 n i 


“ (d)n n! ^ p + c 

n=0 p=() 


V_L_ 

n=0 ^ + ° 


+ 2x 2 exp (—x 2 ) 


d 


iFi ( a-,b-x 2 ) Y 


^(c+1)^ 2 ’ 1 ^ 1 


E 


n=0 ( d+1 )- n! ^ 0 ^ + C+1 


itt (c; d; z 2 ) f; ( ° + ^ x2 " ^-i— 

v ; ( 6 +!)n ra! ^p + a +1 


n=0 


~ ( n \ Jin ™ ^ 1 

+ " F(c+l;d + l;x 2 )]r^^£-L- 

d \tf 0 (b) n n! J^ P + a 


+ ^i-Fi (a + 1; 6 + 1; x 2 ) Y 


(c) n x 2 


n— 1 n 

Evr 


(d) n n! ^p + c 

n=0 7 p=0 


- iF (a; 6; x 2 ) iF (c + 1; d + 1; x 2 ) 

^ iF (a + 1;6 + l;x 2 ) iF(c;d;x 2 ) | 


(60) 


XQ 


While this expression is very difficult to manage by hand, it can easily be processed further 
through a symbolic calculator. 


3.2 Energy-dependent potentials and normalization 

It is well-known that quantum models involving energy-dependent interactions are subject to 
a modified theory [26]. In particular, the usual completeness relation and the corresponding 
F-nornr of wavefunctions associated with energy-dependent potentials change their form. If 
we interpret (JT|) as a Schrodinger equation for an energy-dependent potential V that admits a 
bound state solution u, then the norm N(u) of such a solution must be calculated as follows [26] 

N(u ) = J [1 - V x (t, A)] \u(t, A)| 2 dt. (61) 

D 


Here we use the symbol N rather than the usual notation || ■ ||, because in the strict mathematical 
sense (1611) is not a norm, as it can become negative. Keeping this in mind, we return to the 
Jordan chain dTJ) , ([2]) . According to ([8]), the system is solved by the two functions u and v = u\. 
Let us now calculate the Wronskian of these two functions. To this end, we first find the 
Wronskian’s derivative by taking into account both equations of the Jordan chain, obtaining 


dW, 


u,u x 


dx 


(x, A) = 


= tt(x,A) (u\) xx {x, A) - u x (x, A) u xx (x, A) 

= u(x, A) { — [A — V (x, A)] u\(x, A) — [A — V (x, A)] a u(x, A)} + u\(x, A) [A — V (x, A)] u(x, A) 
= - [1 - V\(x, A)] u 2 (x, A). 





















If we integrate this relation, we obtain the following representation of the Wronskian 

X 

j [1 - V\(t, A)] u 2 (t, A) dt = W Ut u x (x 0 ,X) - W Ut u x (x, A), (62) 

Xo 

for an arbitrary xq € D. We observe that (1621) is a generalization of (1261) when V depends on 
the parameter A. In the next step we make the following observation: if the function u in (1611) is 
real-valued, then we can immediately match the latter relation with the normalization integral 
for systems that feature energy-dependent potentials. Setting D = (xo,x), we get 

N(u) = W U:Ux (x 0 , X) -W UtUk (x,\). (63) 

This identity says that the normalization integral (16 ip can be found entirely through calculating 
derivatives. As mentioned before, this is important in applications, because derivatives are 
generally much easier to find than integrals. Furthermore it is important to point out that the 
identities (1621) and (|63l) neither require the function u to be a bound state solution, nor does 
the associated potential V have to admit such bound states. Both of the latter identities are 
valid independent of any physical meaning, that is, they hold for any solution of ([Tj) and for any 
potential function, as long as the integrals exist. Before we present an example, let us comment 
on the case of a complex-valued function u. Since (f62l) does not contain the absolute value, it is 
not directly applicable. Instead, we can split the integral by 

N(u ) = J [1 — V\(t, A)] Re [u(t, A)] 2 dt + J [1 — V\(t, A)] Im [u(t, A)] 2 dt, 

D D 

note that Re and Im represent the real- and imaginary part of u, respectively. Now, each of 
these two integrals can be identified with our formula (162[) . provided u is replaced by its real- 
and imaginary part. 


3.2.1 Energy-dependent harmonic oscillator system 


Let us now look at a model that features an energy-dependent harmonic oscillator potential. 
Our problem is governed by the following boundary-value problem 

u xx (x, A) + (A — A x 2 ) u(x, A) = 0, (x,A)€M 2 (64) 

lim u(x, A) = 0, A G K. (65) 

laH—Kx> 


This problem is exactly solvable and admits an infinite discrete spectrum (A n ) for an associated 
solution set (u n ), n a nonnegative integer [26], of the form 


— (2 n + 1 Y 


u n (x, A n ) = exp 


-- (2 n + 1) x 2 H n (V2 n + 1 x) , (66) 


where H n stands for the Hermite polynomial of order n [23j. For the sake of simplicity, we will 
restrict ourselves to the simplest case n = 0. We extract from (|66l) 


Aq — 1 


uq(x, A 0 ) 



We will now determine the norm of this function uq by means of the integral (1611) . Afterwards, 
the result is verified by using the differential approach (|63D . The normalization integral ()61l) of 
this function can be evaluated in a straightforward manner 


J [ 1 -(A * 2 )a], a=Ao 


exp (—x 2 ) dx 


J (1 — x 2 ) exp (—x 2 ) dt 
R 
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(67) 


N(uq) 





Let us now demonstrate that the same result can be found using our formula (16311 . First, we 
must determine the partial derivative (uo)\. This involves a little trick, because we are not 
given a solution of ()64l) in terms of A. But since the relation between A (or A n ) and n in (1661) is 
differentiable with differentiable inverse due to nonnegativity of both X n and n, we can interpret 
n as a continuous variable and use the chain rule to obtain 


(u 0 )a(®, A) 


du n 

dn 


(x,X n ) 



|n=0 


du n 

dn 


(x,X n ) 


dXn 

dn 


(n) 


J 171=0 


= -J* 2 «p( 4* 2 )+J exp 


' dH n 

dn 


(a/2 n + 1 x ) 


|n=0 


We incorporate this into the Wronskian in (1631) and arrive after some simplification at the 
following result 


W U0 ,(u 0 )\ 


1 

2 


exp (— x 2 ) [x — H_i(x)\ . 


In the final step we plug this Wronskian into formula (1631) for our norm. Since the numbers 
xo and x correspond to negative and positive infinity, respectively, we must apply limits. This 
yields 


N(uq) 


lim W UiUx (x 0 , A) - lim W UtUx (x, A) 

XQ^X — OO X^XOO 

- lim - exp (-Xq) [x 0 - iL_i(x 0 )] + lim - exp (-x 2 ) [x - -H-i(x)] 

xq — y —oo Z x^xoo Z 

- lim - exp (-Xq) [x 0 - iL_i(x 0 )] 
xq — y —oo z 
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This result coincides with (1671) . as expected. At first sight it seems that the calculations using 
formula (1631) are cumbersome when compared to the simple integration in (1671) . This impression 
is merely due to the simplicity of the present example. If we choose a more general solution of 
(1641) . resolving the integral in (1671) can become very difficult. At the same time, the use of our 
formula (163() will also result in complicated calculations, but it involves only differentiation. 


4 Conclusions 

In this work we constructed a relationship between second-order Jordan chains represented 
in integral and differential forms. This relationship, in a general scenario, allows to find cer¬ 
tain integrals involving functions that are solutions of Schrodinger equations. We applied this 
relationship to obtain conditions to generate regular potentials through the confluent SUSY 
algorithm using the differential representation, to illustrate we used the particle in a box and 
the radial oscillator systems and in both cases it was shown how to find integrals related to the 
involved special functions of each system entirely through derivatives. It is worth mention that 
the exactly solvable SUSY partner potentials of the radial oscillator presented in this work as 
example had not been reported. Also, since the limits of the found integrals are not fixed, they 
can be used to calculate probabilities in any arbitrary interval inside the domain of definition 
of the potential. Moreover, the relationship was used to find normalization constants of wave 














functions of energy dependent potentials in quantum mechanics, it was exemplified with the 
ground state of an energy dependent harmonic oscillator. The examples explored in this article 
are far from being exhaustive but have rather exemplary character since for every quantum 
system integrals with different special functions can be analyzed 
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